Post-Einsteinian tests of linearized gravitation 
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The general relativistic treatment of gravitation can be extended by preserving the geometrical 
nature of the theory but modifying the form of the coupling between curvature and stress tensors. 
The gravitation constant is thus replaced by two running coupling constants which depend on scale 
and differ in the sectors of traceless and traced tensors. When calculated in the solar system in a 
linearized approximation, the metric is described by two gravitation potentials. This extends the 
parametrized post-Newtonian (PPN) phenomenological framework while allowing one to preserve 
compatibility with gravity tests performed in the solar system. Consequences of this extension are 
drawn here for phenomena correctly treated in the linear approximation. We obtain a Pioneer-like 
anomaly for probes with an eccentric motion as well as a range dependence of Eddington parameter 
7 to be seen in light deflection experiments. 

PACS: 04.20.-q, 04.80.Cc 



I. INTRODUCTION 

Experimental tests of gravity performed in the solar 
system show a good agreement with General Relativity. 
This statement can be put under a quantitative form by 
usingthe parametrized post-Newtonian (PPN) formal- 
ism [l| or by bounding deviations of the gravity force 
law from its standard form General Relativity is 
however challenged by observations performed at galac- 
tic or cosmological scales. Due to the good agreement 
of tests with General Relativity, the anomalies seen in 
the rotation curves of galaxies or in the relation between 
redshifts and luminosities for supernovae are commonly 
accounted for by introducing dark matter and dark en- 
ergy components designed to this purpose 0, 0, • As 
long as these new components are not detected by other 
means, the anomalies can also be thought of as conse- 
uences of modifications of gravity laws at large scales 

An important requirement to be met by such modifi- 
cations is that they remain compatible with the bounds 
set by solar system tests. The recently observed Pioneer 
anomaly might be a key piece of information in this con- 
text by pointing at some anomalous behaviour of gravity 
at a scale of the order of the size of the solar system [ij . 
The effect appears as an anomalous acceleration recorded 
on Pioneer 10/11 probes during their flight to the outer 
solar system. The measured Doppler tracking data are 
compared with predictions of General Relativity and the 
residuals then expressed as an unmodelled acceleration 
for the probes. This anomalous acceleration is directed 
towards the Sun with an approximately constant ampli- 
tude 8 x 10 -10 ms -2 over a large range, 20 to 70 astro- 
nomical units (UA), of heliocentric distances |10| . 
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Up to now, the anomaly has escaped all attempts of 
explanation as a systematic effect generated by the space- 
craft itself or its environment, and it has no more been 
convincingly derived as a consequence of new physics (see 
the references in ^3,^3)- This status should motivate 
further scrutiny of any potential origin of the anomaly. 
The aim of a better understanding of deep space nav- 
igation is by itself a sufficient motivation and the 
possibility that the Pioneer anomaly be the first hint of 
a modification of gravity laws at large scales cannot be 
let aside investigations |14|. 

The long standing opposition between General Relativ- 
ity and Quantum Theory has led to the view that the for- 
mer can only be an effective theory of gravity which could 
be altered at short or large length scales. Radiative cor- 
rections naturally lead to an immersion of General Rela- 
tivity within a class of fourth order theories 0, 0, 
with potential effects at short ranges Modifications 
may also be expected to appear at larger length scales 
[ill l20j with implications in astrophysics and cos- 
mology [23, 123I |2J| . In contrast to Einstein theory which 
is natively non renormalizable |25j , fourth order theories 
show renormalizability as well as asymptotical freedom 
at high energies [pej [2(j. The extension of gravitation 
theory at scales not already constrained by experiments 
can thus be tackled by studying renormalization group 
trajectories p7j . 

Renormalizability of the family of fourth order theories 
however comes with a counterpart, namely the problem 
of unitarity associated with ghosts which is often thought 
to flaw fundamental field theories. It has however been 
convincingly argued that it does not constitute a defini- 
tive deadend for an effective field theory valid in a lim- 
ited scale domain psj . In particular, the departure from 
unitarity is expected to be negligible at ordinary scales 
tested in present day universe |2g, leading to sensible 
calculations for observable phenomena |3fl l3ll l32l ] . 

In the present paper, we investigate the potential con- 
sequences of modifying the gravitation theory at length 
scales of the order of the size of the solar system [33j . 
We show that there exist natural extensions of General 
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Relativity which preserve its essential geometric features 
such as the identification of gravitation fields with the 
space-time metric, the definition of geodetic motions and 
the equivalence principle. Only the form of the coupling 
between the Einstein curvature tensor and the energy- 
momentum tensor is modified. In agreement with be- 
haviours observed for radiative corrections of General 
Relativity, we assume that this coupling can become scale 
dependent and differ in the two sectors of traceless and 
traced components. As the application of the new frame- 
work is primarily focused on potential consequences in 
the outer solar system, we then treat the gravitation the- 
ory in a linearized approximation and describe the grav- 
ity source as punctual, isotropic and static. 

We thus show that the metric is characterized by two 
potentials directly related to the two running coupling 
constants which take the place of the single Newtonian 
gravitation constant. The first potential generalizes 
the usual Newton gravitation potential and its deviation 
from the latter suffers stringent bounds deduced from the 
remarkable agreement of General Relativity with gravity 
tests. Meanwhile, the second potential opens addi- 
tional freedom offering the possibility to accomodate a 
Pioneer-like anomaly for probes with eccentric trajecto- 
ries. The second potential $p may also be thought of as 
promoting the Eddington parameter 7 to the status of 
a range dependent function. It thus leads to the predic- 
tion of other anomalies to be recorded on null geodesies, 
in particular in time delay, deflection or Doppler track- 
ing experiments performed on probes passing behind the 
Sun [MS HI III HE! • These anomalies which have 
the same origin as the Pioneer-like anomaly in the new 
framework will be discussed in detail below. Hopefully, 
their confrontation to observations will either provide us 
with data points to be compared with the Pioneer data 
points or lead to new bounds on deviations of alternative 
gravitation theories from General Relativity. Note that 
the effects associated with non linearities of gravitation, 
that is also with the Eddington parameter /?, will only 
be discussed briefly since they need not be considered in 
situations studied here (more details on this point below). 

The following sections begin by a description of the 
proposed extension of Einstein gravitation theory in a lin- 
earized approximation. Considering a punctual, isotropic 
and static gravitational source, we deduce the expression 
of the modified metric in terms of two potentials <I> n and 
$p. We then examine the consequences of the modifi- 
cations on the third Kepler law, on Doppler tracking of 
probes and on light deflection. We finally use the sim- 
plest example, which involves only three additional pa- 
rameters, as a benchmark for comparing the predictions 
of the modified framework to available observations. 



II. MODIFIED GRAVITATION EQUATIONS 

We first introduce the extension of Einstein theory 
by discussing its motivations and writing down modified 



post-Einsteinian gravitation equations. The term 'post- 
Einsteinian' has a twofold meaning alluded to in the In- 
troduction. It first means that the geometrical nature 
of gravity, the very core of Einstein theory and, further- 
more, one of the best ever tested properties of the physi- 
cal world, is left untouched. But it also points to the fact 
that the coupling between curvature and stress tensors 
suffers a modification from its Einsteinian form. 

As a first step in this discussion, we emphasize that 
the geometric foundations of the theory are preserved : 
the gravitational field is identified with a metric tensor 
g^ Vl the motions are described by associated geodesies 
and, as a consequence, the equivalence principle is left 
untouched. As is well known, this principle is verified at 
distances ranging from the millimeter in laboratory ex- 
periments ( 4(| and references in) to the sizes of Earth- 
Moon orbit ^3 or Sun-Mars orbit 0, • Tne relative 
accuracy of these tests, better than I0~ 12 for some of 
them, is good enough to discard any interpretation of the 
Pioneer anomaly from a violation of the equivalence prin- 
ciple ^(j . This statement does not entail that the equiv- 
alence principle is exact, but it means that the potential 
deviations from General Relativity studied here have a 
larger amplitude than the violations of the equivalence 
principle expected from unification models |44l [45j . This 
is why they can be analyzed in a metric theory sharing 
the conceptual basis of General Relativity. 

As already discussed, we focus the attention on the lin- 
earized theory where the gravitation field is represented 
as a small perturbation h^ v of the Minkowski metric r\p V 

9nv = V^u + (I) 
rif,, v = diag(l, -1, -1, -1) , <C 1 

The field h^ v is written as a function of position x in 
spacetime or, equivalently, of wavevector k 



''"" U) 1 J^ e ~ ikXh »» [k] 



(2) 



Riemann, Ricci and scalar curvatures are easily written 
in momentum representation, at first order in h^, 



R 



k\k u h^p kxkph^u k^k v h\p + kpkph\ v 



k^h^p kpk a h a p kpk a h a ^ + k^kpTj^h^ 



R — k tj^ 1 h^i/ k^k hp^ii 



(3) 



We use the same sign conventions as in [4(j , indices being 
raised or lowered using Minkowski metric Tjp U . 

These expressions can be recovered from Einstein cur- 
vature tensor 



E pv — Rpv ^^jfiuR 

R^v — Efj, v — -rjfivE , E = —R 



(4) 



R\l 



k\k u R^p k\kpR^ LV kpk v R\p + kpkpR\ v 
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Ep V may be written in terms of transverse projectors 



E„ 



(5) 



T^lxv — Tlliv ^2 ' ^ — k p kp , k p TX^ v — 

It may then be decomposed as the sum of two indepen- 
dent components with different conformal weights 

Epv = E + E 1 



E° = 7T° E pa = -7T° k 2 h pa 
^ixu — n jj,vpa- LJ 2 w 

^fiv — " fiiipcr-'-' "p.iypcr^' " 



(6) 



where 7r° and tt 1 project transverse tensors onto their 
traceless and traced components 



^Lip^va \ ^iia^vp i 

— 7T 

2 fJbUpu 

7^ pyTT p a 



pvpa g 
° TT Pa = 

pVpO 



(7) 



As discussed in the introduction, we disregard the ef- 
fects of rotation and non sphericity of the Sun which 
affect significantly the motions of inner planets but have 
small effects on motions in the outer solar system or on 
propagation of light. More precisely, we consider that 
these effects are properly taken into account in the stan- 
dard description based on General Relativity (see for ex- 
ample the calculation of Doppler tracking data for Pi- 
oneer probes in |lfj|) and that they have the same im- 
pact in the modified and standard descriptions. As a 
consequence, the anomalies evaluated by subtracting the 
standard result from the modified one can be calculated 
with the assumptions of stationarity and isotropy. The 
same statements can be applied essentially unmodified 
to the assumption of linearity (see more details on this 
discussion in the next section). 

Fields are written as functions of spherical coordinates 

x p = (ci, r cos # cos i^9,r cos # sin iy9,r sin #) (8) 

where c denotes light velocity, t and r the time and 
radius, 9 and if the colatitude and azimuth angles. 
Isotropic quantities are functions of x° and r only. Then, 
the Einstein tensor © is read in terms of two potentials 



Ep U = TT 



/ii/OO 



2^0 



/if 00 



2fc 2 $ 1 



(9) 



When stationarity is also assumed, these relations are 
simplified due to the fact that ko = and k 2 = — k 2 . 

The Einstein tensor E^ v is transverse as a consequence 
of Bianchi identity and this is also the case for the stress 
tensor Tpy, due to energy- momentum conservation. The 
second condition also implies geodesic motion for test 
masses, in consistency with the geometrical interpreta- 
tion of gravity which is preserved by the extension of 
Einsteinian theory studied here. In General Relativity, 
these two tensors are merely proportional to each other 

M 



Efj, v — 8ttGnTp V 



(11) 



with Gn the Newton constant. Now this equation may 
easily be generalized while preserving transversality so 
that, in the linearized theory, the relation between the 
two tensors takes the form 



Et. 



XpvXp 



rp\p 



a = 0,1 



(12) 



The response functions Xuv are constrained by the 
transversality condition but they may depend on momen- 
tum and differ in the two sectors, thus extending the form 
of the coupling assumed in Einsteinian theory. 

We treat the Sun as a static gravitational mass M 
localized at the origin x — of the coordinate system 



T M „[k] = r]p r)„oMc 2 



(13) 



We have introduced a notation for the Fourier transform 
of stationary quantities 



Tp4k]=Tp4k}27r6(k ) 



(14) 



The potentials thus have stationary and isotropic expres- 
sions given by Poisson like equations 



k 2 $ a [k] = 



2G a [k]M 



a = 0,1 



(15) 



With the assumptions of stationarity and isotropy, and 
for a pressureless gravitational source (see ea. (ll3ll '). the 
two running constants G a constitute the information on 
the linear response functions x a which can be extracted 
from experiments 



G°[k] 

c 4 
&\k] 



gXoooo[0,k] 



-X 



/ji/00 



[0,k] 



(16) 



$° and <I >1 are functions of x° and r or, alternatively 
in momentum representation, of ko and |k| with k 2 — 
fcp — k 2 . Conversely, $° and "5 1 can be written in terms 
of the traceless and traced curvatures © 



2$ °-^(kV^° 



2$ x 



k 2 1 
3 (k 2 ) 2 



mi 



(10) 



Einstein theory corresponds to the simple case where the 
two running couplings reduce to a single constant Gn- 
Note also that these two running constants G a do not 
exhaust the information contained in the linear response 
functions \ a ■ It is f° r example easy to show that equation 
(|12f> allows for three independent scalar functions while 
still preserving transversality, stationarity and isotropy. 
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But the third one cannot be tested in the solar system 
due to the absence of pressure terms in the gravitational 
source l|T3jl . 

In order to underline the potential deviations from 
General Relativity, we will introduce a generic notation 
showing the separation between a standard expression 
[ ] s t (* e ^ ne prediction of General Relativity) and an 
anomaly indicated by the symbol S. For example, the 
two potentials are the sums of a standard expression 



[*V)L = [*V)L = 



G N M 



') (17) 



and of a deviation which has to remain small in order to 
fit existing gravity tests 

S°(r) = [$ a (r)] st + 5$ a (r) , |<5$ a (r)| < 1 (18) 

The symbol <C is used here to suggest that the quantities 
6& a are smaller than the Newtonian potential <f> which 
is already smaller than unity in the linearized approach. 
Since these functions may have different dependences ver- 
sus r, this can only be a qualitative statement at the 
moment. The main objective of the present paper is to 
discuss the potential phenomenological consequences of 
<5<f> a in the range of momentum or, equivalently, of length 
scales tested in the solar system. 

With this purpose in mind, we will avoid favoring any 
particular functional dependence of the anomalous run- 
ning constants SG a or of the associated anomalous po- 
tentials SQ a . In order to illustrate the discussion, we 
will however use occasionally as a "benchmark" a spe- 
cific form of the potentials 



$ Q (r) 



G a M ( a Mr 



a = 0, 1 



(19) 



This form superimposes linear terms proportional to r 
to the standard terms proportional to 1/r. Equivalently, 
the running constants are the sums of standard constant 
terms and of infrared corrections proportional to 1 /k 2 



G a [k] ~ G a 



2( a 



0,1 



(20) 



The main advantage of this simple example is that the 
phenomenology is now determined by four constants, the 
interpretation of which is given below. In fact constants 
of integration appear when solving equations i|15|) with 
the running constants Ij20|l . They give rise to terms pro- 
portional to r° and r 2 besides those written in (|19l) . Since 
constant corrections to the metric do not affect curva- 
tures while quadratic terms correspond to constant cur- 
vatures relevant at larger galactic or cosmological scales, 
both can be ignored in the solar system. 

The linear dependence of the terms proportional to 
C" in (|19fl also calls for remarks. Should it hold over 
galactic length scales, this dependence would induce ef- 
fects incompatible with observations. This means that 
(|19fl can only be considered as effective potentials valid 



in restricted ranges of distances or momenta. The lin- 
ear terms can for example result from an expansion at 
rCrA of Yukawa type potentials (1/r) exp(r/r\) which 
are regular at large distances 48] . Note that radia- 
tive corrections to General Relativity are described by 
modifications of the lagrangian which behave, at low- 
est order, as quadratic forms of the curvature tensors 

El 13 El E3 

and then give rise, in the linear approxi- 
mation, to Yukawa corrections. Though it is not possible 
to meet the goal of accomodating the Pioneer anomaly 
in a frame compatible with planetary tests by ad ding 
Yukawa corrections to the single Newton potential |48| . 
the presence of two potentials with different values in the 
two sectors now opens an additional freedom. 

In the following, we consider modifications of the grav- 
ity theory described by two potentials <I> a (r). Occasion- 
ally, we use the simple forms (|19ll or (|20|l for deriving 
preliminary predictions for experiments in the solar sys- 
tem, keeping in mind that a large distance regulator is 
implicit when discussing for example the deflection of 
light coming from stars or extragalactic sources. 



III. MODIFIED METRIC TENSORS 

We now write the metric obtained as a solution of the 
modified equations (|15I16(I and to be used for discussing 
gravity tests in the forthcoming sections. 

To this aim, we adopt the gauge convention of the PPN 
formalism lj with the metric written as a stationary and 
isotropic function of spherical coordinates © 



ds 2 = go c 2 dt 2 + g rr (dr 2 + r 2 (d6 2 + : 



l )) (21) 



This gauge choice fixes the expressions of the metric com- 
ponents in terms of the two potentials i|10|) 



.goo = 1 + 2$at , $at 



4$° _ (J) 1 



(22) 



2($° - $!) 

g rr = -1 + 2$ N -2$ P , $ P = ^ '- 

$jv and <I>p have been introduced as two linear combi- 
nations of <E>° and $ 1 . <E>at is defined from the difference 
(,<?oo — 1) and thus identified as a generalization of New- 
ton potential 4> while <f>p is defined from (— gooSrr — 1) 
and measures the difference between the two sectors of 
traceless and traced curvatures. 

As $° and $ 1 , and <&p obey Poisson equations with 
running constants Gn and Gp written as linear combi- 
nations of G° and G 1 



k 2 $ a [k] = 



2G a [k]M 



- _ 4G° - G 1 ~ 

= „ ; LrP 



a = N,P (23) 
_ 2(G° - G 1 ) 



3 ' 3 

Einstein theory corresponds to the standard expressions 

G N M 



[*Jv(r)L = #r) 



, [$p(r)] st = (24) 
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and deviations have to remain small 

$„(r) = [*a(r)] rt + <y« (r) , |**a(r)| « 1 (25) 

According to the discussion of the preceding section, 
we will occasionally use the specific form (|19|l of the po- 
tentials 

G a M CMr , s 

with the new cocfficnts obtained as in (|2U[I . With the 
simple model (I26|) . the phenomenology is determined by 
four constants, the Newton constant Gpj and three small 
parameters Gp, (n, Cp- This allows one to reach more 
specific conclusions with however the drawback of a loss 
of generality with respect to i|25|) . This is why we will 
as far as possible use the general form (1221) of the metric 
and restrict the use of the simple model (|26[) to a few 
specific occasions. 

At this point, it is worth comparing the metric 122|) 
with the parametrized post-Newtonian (PPN) expression 
introduced by Eddington and then developed by sev- 
eral physicists [H Ull El El . The Eddin g ton PPN met- 
ric is read in the isotropic gauge lj2"Tl) 

.900 = 1 + 2a0 + 2/3(j? + ... 

g rr = -l + 2 7 cj)+ ... (27) 

As the three Eddington parameters are constant, the first 
one a is usually eliminated through a redefinition of the 
Newton constant. We show now that the metric 1221) can 
be regarded as similar to l|27|) but with the quantities a, 
(3 and 7 promoted from the status of a constant to that 
of a function. To this aim, we rewrite l|27(l as the sum of 
a standard expression obtained for a = = 7 = 1 and 
of a deviation 5g^ v 

5 5oo ~2(a-l)0 + 2(/3-l)0 2 + ... 

Sg rr ~ 2(7 - \)<j) + . . . (28) 

and we identify the terms appearing in (|28|l with the 
anomalous potentials in l|25|) . 

First, the term (en — l)(f> in l|28[l becomes the function 
<5<f>Ar. It may have a more general r-dependence which 
also means that a is no longer a constant. Similarly, the 
potential $p is another function of r which somewhat 
generalizes the term (7 — 1)0 in (|2"81) . As a consequence, 
gravity tests expressed as an upper bound on the con- 
stant (7 — 1) in the PPN formalism will become a bound 
on the function $p in the new framework. In loose words, 
the bound on the Eddington parameter 7 may now be- 
come non uniform at different distances in the extended 
framework. This statement will be made more quanti- 
tative in the following. We will see that a long range 
variation of $p can lead to a gravitational interpretation 
of the Pioneer anomaly while circumventing the main ob- 
jection to such an interpretation, namely the absence of 
a similar effect on planets ^(| . The presence of the sec- 
ond potential $p at short heliocentric distances will also 



affect time delay or light deflection experiments, opening 
the possibilit y to see it by reanalyzing data of deflection 
experiments |38| . 

The manifestation of the second potential $p in 
anomalies usually ascribed to (7 — 1) is our primary sub- 
ject of interest here. Before focusing the discussion on 
this subject, let us open a parenthesis on the non lineari- 
ties of gravitation theory described by the parameter p in 
the PPN metric l|28|l . These effects are known to be sig- 
nificant in the inner part of the solar system by affecting 
the perihelion precession of Mercury [54j as well as the 
polarization of the orbit of Moon around Earth induced 
by the Sun [5^] . In contrast, they are small for probes 
in the outer part of the solar system or for propagation 
of light and this is a first good reason for disregarding 
them in the context of the present paper. Let us briefly 
sketch how this simple argument can be put in a more 
quantitative form. 

The expansion 1221) of goo can be pushed one order 
further to obtain an expression comparable to the PPN 
expression (|28|1 

Sgoo = 2*$$ + 25$% ] (29) 

is just the first order potential discussed up to now 
(2) 

and & N the second order contribution to the metric com- 

(2) (2) 
poncnt goo- Then S& y N ' is the anomalous part of & y N ' 

which takes the place occupied by (/3 — l)cf> 2 in the PPN 
expression (|28() . In loose words, 6&ff promotes a con- 
stant anomaly (/3 — 1) to the status of an extra function. 
This function should play a significant role in the evalu- 
ation of the Mercury perihelion precession or of the orbit 
of Moon, a discussion which requires a non linear treat- 
ment of gravitational perturbations and is not done in 
the present paper. However we want to take care of its 
existence if only to make clear that it does not spoil the 
conclusions of the linear treatment. 

To this aim, we now consider an argument which plays 
a key role in the presentation of the PPN metric l|28|l . 
given up to first order for g rr but to second order for 
<7oo- The idea is that terms proportional to (3 may have 
the same impact on motions than those proportional to 
7, though they appear at the next order in the metric. 
The spatial metric components are indeed multiplied by 
the kinetic factor v 2 /c 2 in the action of test masses and 
this factor has the same magnitude as \4>\ for virialized 
motions, in particular for circular motions. Now, as a 
consequence of the same argument, the effects associated 
to 7 should dominate those associated to (3 as soon as 
the kinetic factor v 2 /c 2 is much larger than \<p\. This 
is the case for the two situations studied in the present 
paper, namely the motion of Pioneer-like probes escaping 
the solar system and the deflection of light. In other 
words, this argument allows one to delineate the category 
of phenomena correctly accounted for in the linearized 
treatment and the problems studied in this paper are 
precisely chosen within this category. 
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IV. MODIFIED GEODESICS 

The geometric foundations of General Relativity are 
left unchanged in the new framework and motions are 
described by geodesies associated with the stationary and 
isotropic metric Q22J1 . We characterize these geodesies 
and, as a first application, we evaluate the modification 
of the third Kepler law. 

The motion of a probe may be written down as the 
Hamilton- Jacobi equation for the action S, that is equiv- 
alently the phase of the associated field, 



Another form of the equation of motion is the geodetic 
equation obtained by differentiating H30|) or (|34[) 



Du^ du 

~ds~ ~ ' " 
1 



M - -IT ~ ^^ P u v vP = 



(36) 



{dv9 w + d p g^ - d^g vp 



D is the covariant derivative and r Mj „ p the Christof- 
fel symbols associated with the metric. This equation 
can also be written in terms of accelerations defined as 
Christoffel symbols projected along the motion 



g^daSduS = m 2 c 2 



(30) 



g^ v is the inverse of the metric and m the mass of the 
probe, with m = for light. This equation is easily 
solved for a stationary and isotropic metric |46j ] with no 
dependence on time t and angle variables ip and 9 

dtg^v = d^g^ = deg^ = (31) 

Energy E and angular momenta J and K are conserved 

dt 



E 



J 



K 



-d t S = mc g 00 
-d v S = mcg vip 
-d e S = mcgee 



ds 
dip 
ds 
M 

ds 



(32) 



and the action S can be written in terms of a phase shift 
function S r depending only on r 



S = -Et + Jip + S r 

dfS — dipS^ — - yricQw - -! — 
as 



(33) 



Without loss of generality, the trajectory has been as- 
sumed to take place in the plane 6 = ? (K = 0). 

The Hamilton-Jacobi equations H3()|l also expresses the 
normalization of the 4-velocity u M 



u; 
E 2 



dx v 
ds 

9oo [ 1 



g^u^u" = 1 

1 J 2 

g rr m 2 c 2 r 2 



(34) 



9oog r 



With the metric (|22|l . this relation is read as follows in 
terms of the potentials 



E 2 J 2 
-5-j = l + 2$ Ar + (l + 4ci> )- 5 - 5 - 



-" (J 



(35) 



dUf, 
ds 



p = r v » v p 



(37) 



In particular, the radial acceleration can be rewritten 



(1 - 2$ 



d$ N d$ P / dr 



dr dr \ ds 



2 dr \ m z c z i 



(38) 



In the following, we use this expression linearized at first 
order in the gravitational perturbation and written as the 
sum of the standard expression and of an anomaly 



IV,ss 

[r r ,ss] st 



Fry 



dr 



20 -2r-f 



J 2 



d5$ 



N 



dr 



1 + 2(5$at - 2r 



dr J m 2 c 2 r 3 
dS^ N 



dr 



d (it. J2 



d$ P f dr 



dr \ ds 



(39) 



This expression makes it easy to identify the contribu- 
tions of (5$tv and $p to potential anomalies to be seen 
in the weak field approximation. The role played by the 
contribution of 4>p in Pioneer like anomaly and deflection 
tests is discussed in detail in the next sections. 

As a first application, we now evaluate the modification 
of the third Kepler law for circular orbits which depends 
essentially on 5$jv- To this purpose, it is convenient to 
introduce the squares of the radial and angular velocities 
with respect to time, v 2 and v 2 respectively, 



2 _ -2 

v r = -g rr r 



2 2-2 
« w = ~9rrT If 



_ dr 
= ~di 

. _ dip 

ip "~dt 



(40) 



We now differentiate the normalization relation (|34|l with 
the constraints that E and J are conserved and obtain a 
generalized Kepler law 



Note that (l + 4$°) is the linearized form of the confor- 
mally invariant quantity {—goo/ 9rr) (see ea. l|22|) '). 



1 



2 (1 - \rd r g rr ) 



d r goo + glodr 



( -I 



\9ooC 2 , 



(41) 
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For circular orbits 



0), this law is read 



from the emitter to the receiver 



rd r g 0Q 



(42) 



For weak gravitational fields and at leading order, this 
Kepler law depends only on g 00 , that is also on <&jv, 



2 • 2 

r ip 



c 2 rd r $ 



N 



(43) 



An anomaly 5&n of leads to an anomaly in the third 
Kepler law for a circular orbit, as is made evident by 
decomposing in standard and anomalous parts 



KL+K 



VI st 



fist 

z 2 rd r 



(44) 



5vl ~ c 2 rd r d^> 



N 



The agreement between General Relativity and tests per- 
formed on planets tells us that d r 8$N has to be small. In 
any case, it is certainly too small to explain the Pioneer 
anomaly (see the discussions in and 0])- Conse- 
quently, the deviations from General Relativity studied 
in the present paper are mostly ascribable to the effects 
of the second potential <i>p discussed in the next sections. 

Note however that more accurate statements taking 
into account higher order effects can be obtained by ap- 
plying exact expressions written in the present section to 
the metric. According to the argument presented at the 
end of section lrnl this analysis should let the conclusions 
of the present paper unaffected. 



V. DOPPLER TRACKING 

We show now that the second potential <£>p, which does 
not affect appreciably the third Kepler law for circular 
motions, may in contrast lead to a Pioneer like anomaly 
for probes having an eccentric motion. To this aim, we 
perform within the extended framework the calculation 
of Doppler tracking of such probes. We do not enter 
into the technical details described in the analysis of the 
Pioneer anomaly in |lf)j |. 

As a first step, we write the expression of the Doppler 
frequency shift and its time variation in terms of the 
two gravitation potentials and of the probe velocity. The 
tracking is built up on lightlike signal propagating from 
Earth to the remote probe and back from the probe to 
Earth. Earth follows a nearly circular geodesic whereas 
the probe follows an eccentric geodesic escaping the so- 
lar system. The up-link radio signal leaves the Earth 
at the location x7 = (t7 , ) and attains the probe at 
%2 = (^2,X2). There, it is transponded into a down- link 
radio signal leaving the probe at X2 and reaching the 
Earth at xf = (t^,x+). The two links are null geodesies 
with their endpoints, corresponding to emission and re- 
ception, on the geodesies of the station on Earth and the 
probe. They are described by a function T which mea- 
sures the time taken by the lightlike signal to propagate 



n T T(X! ,X 2 ) 

i 2 +T(x+,x 2 ) 



(45) 



This function T is calculated by considering the up and 
down links as the advanced and retarded parts of a light- 
cone originating from a; 2 . As there is only one lightlikc 
path from a spatial point to another in the weak gravi- 
tational field of the solar system, this definition is unam- 
biguous. For the metric 112211 furthermore, the function T 
is symmetric in its two arguments. Its explicit expression 
is derived from the action (|33|) in appendix 1X1 

The endpoints xf and X2 follow geodesies with veloci- 
ties uf 1 and U2 ■ Since they are related by equations 145(1 , 
the corresponding proper times ds7 , ds^ and ds2 obey 

(u7 ■ n7) , , (uf ■ nf) , , , , 
1 -dsT = ds 2 = — y-dst (46) 



("2 • n 2 ) 



(u 2 ■ n+) 



ri^ represent the directions of wavevectors at endpoints 
a = 1, 2 and ■ denotes a 4-dimensional scalar product 



( n t)o = -(nf)o = 1 
(u ■ n) = u^n^ 



(47) 



The up-link signal leaves the Earth with a frequency 
determined by comparison with a clock located at the 
emission station, it travels to the probe where it is 
transponded and multiplied by a constant factor. It trav- 
els back to the Earth where the frequency of the received 
down-link signal is compared with a clock located at the 
reception station. The Doppler observable is defined from 
the ratio of the emitted and received frequencies after 
correcting for the known constant factor. The frequen- 
cies are measured against clocks which tick according to 
proper times ds7 and dsf , so that the round trip Doppler 
shift y is built up on the ratio of these quantities 



ds 1 



(48) 



Although only the linear part will be needed, this precise 
convention has been chosen for the sake of simplifying 
forthcoming expressions. As ds7 and dsf obey the 
round trip Doppler shift is finally given by 



{u 1 ■ n 1 ) (u 2 ■ TI2) 
(u 2 ■ rq) [ut ■ n+) 



(49) 



This is a properly defined observable which is explicitly 
invariant under gauge changes. It must also be empha- 
sized that it accounts for the effects of motions (Doppler 
effect) and of gravitation (Einstein redshift effect). 

The Pioneer observable is then obtained by looking at 
the Doppler shift (|49() as a function of time. Following 
|10| . we convene to write it as an acceleration a with 



2ads 



-c dy 



(50) 
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At this point, ds represents the proper time delivered by 
a reference clock to be specified later on. Note that the 
definition of a thus depends on the choice of this refer- 
ence clock whereas the definition of ads does not. We 
also assume the observers and probe to follow geodesies, 
that is motions with vanishing covariant acceleration 
Du\ = Duf = Du 2 = (corrections associated with 
the position of stations on Earth are taken into account 
separately). Using the expression 149(1 . we thus obtain 
the relation 



2ads (u 1 ■ Dn 1 ) (itj~ • Dnf) 

c 2 (%•%) (uf-nt) 

{u 2 ■ Dn 2 ) ^ (u 2 ■ Dn^) 
(u 2 ■ n+) (u 2 ■ n 2 ) 



(51) 



This relation is exact in the context of our simplifying as- 
sumptions. As already emphasized, the small corrections 
associated with perturbations not taken into account in 
this context are supposed to be linearly superposed to the 
anomaly and, therefore, to disappear in the forthcoming 
expressions obtained after a subtraction. 

The Doppler acceleration lj5T|l has a twofold depen- 
dence versus the metric. First, the Christoffcl symbols 
which enter the covariant derivatives have their devi- 
ation from standard form given by (|39(l . Second, the 
time delays which enter the expressions of the anoma- 
lous wavevectors are deduced from (I47|l as 



(*n±) = 



0(xf 



(52) 



ST is the anomalous part of the time delay function and 
is computed in [X] One deduces the variations of the 
Doppler shift lHHJ 



ds+ 
dsi 



(u 2 - [n 2 ] st ) (uf- [n£]J 



Sy = - 



I st' V 1 11 J St' 

1 ccMT(x^,x 2 ) 



(u 2 



I si 



ds 2 

cd<5T(x+,x 2 ) 



(u 2 



dso 



(53) 



time delay and Christoffel symbols respectively 



5 a 

2SaTds 2 
25a^ds 2 



Sap + <5ap 
-dSy 



(54) 



ds 



ds 2 



1 K - n i 

(<$r S2S2 



ds 



(u 2 ■ r4) 



(uj ■ nj) 
{5T S2S2 ■ rq) 
(u 2 ■ n%) 



We have now specified the reference clock to be on the 
probe (ds = ds 2 ) but another choice would not affect this 
expression at the first order in the gravitation potentials. 
We have also replaced all quantities but ST and ST by 
their zeroth order approximation. 

Using equations ((38(1 , one sees that the contributions 
of the Christoffel symbols almost compensate in (154(1 for 
the Earth which has a nearly vanishing radial velocity. 
Contributions from the Earth motion may then be ne- 
glected so that the anomalous acceleration 1(54(1 may be 
rewritten, at leading order, 



25a 



P 



25a r 



ds\ 
-25T r 



c5T(xr,x 2 ) + C (5T(x+ x 2 )) 



(55) 



Inserting in equations ((55|) explicit expressions of time 
delays l(A10() and Christoffel symbols 1(391) while neglect- 
ing terms proportional to the angular velocity of the 
probe, one obtains the two contributions to the anomaly 
in terms of the anomalous potentials 5$>n and <5<f>p 



Sap 
Sa-p 



*((28*„-6* P )p. 
ds 2 \ ds 2 



d5$ N d5$ P ( dr 2 
dr 2 dr 2 \ ds 2 



(56) 



It was shown in the preceding section that the anoma- 
lous part 5<&n of the first potential <&at had to remain 
very small to fit planetary data. This implies that the 
potential anomalous acceleration is mainly ascribed to 
the second potential <E>p. We will see in the section IVlII 
that the terms 5ar and Sar thus have similar contribu- 
tions. Incidentally this will prove that none of the two 
contributions could be disregarded. 



VI. LIGHT DEFLECTION 



Collecting these results, we obtain the anomalous ac- 
celeration 5a by subtracting the value calculated in stan- 
dard Einstein theory (<5$p = 5$n = 0) from the mod- 
ified value calculated with 5$p ^ or 5$>n ^ 0. This 
quantity 5a corresponds to the Pioneer anomaly denoted 
by cip in ^} • We write it as the sum of two contributions 
Sap and Sar corresponding to the anomalies induced by 



As argued in section IIIII the linearized treatment is 
sufficient for probes having a large kinetic energy and 
this certainly includes the study of light waves. We now 
study the effect of the second potential $p on Eddington- 
like light deflection experiments 0, building up the de- 
scription on the treatment of the time delay function and 
Doppler shifts presented in section Ivl 
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There are indeed close similarities between Doppler 
tracking of spacecrafts and deflection experiments and 
the most recent deflection experiments use Doppler track- 
ing techniques [3^, |5(j . As an important difference how- 
ever, the effect of the time delay is favored by a kind 
of amplification during occultations by the gravitational 
source which play a key role in deflection experiments. As 
a consequence, the modifications of motions of endpoints 
are found to play a negligible role in the observables. We 
show below that, in these conditions, the deflection of 
light is described by an effective Eddington 7 parameter, 
which is now a function rather than a constant. The usual 
expression of deflection can be used with the function 7 
evaluated at the impact parameter p, that is the closest 
distance of approach of the light ray to the gravitational 
source. 

We consider a Doppler measurement on a down-link 
radio signal from the probe to a station on Earth. We 
use the same notation as in the preceding section with 1*2 
and Ui denoting the space-time velocities of the emitter 
and receiver, n\ and n\ the directions of wavevectors 
at emission and reception, ds2 and ds^ the respective 
proper times. We define the Doppler shift y+ on this 
one-way tracking technique by using (|46|l 



ds 2 



{u 2 



(ut-nt) 



(57) 



As previously, this formula describes Doppler as well as 
Einstein effects. For simplicity, we suppose the emitter 
to be at rest at a radius T2 larger than the Sun-Earth 
distance r\ — iAU (with angle (^2) and the receiver on 
Earth to follow a circular orbit with constant radius n 
(with angle ipi). The angular velocity ip± of Earth is given 
by the third Kepler law (see section I1V|) . The Doppler 
shift is then deduced from the time delay function T 



e~ y+ = K 



K = 



dT 
~d~h 



(58) 



y/gooin) + g rr (ri)ri 2 tpl/c 2 

K is the ratio of rates for the emitter and receiver clocks. 
This factor was unity in the two-way technique studied 
in the preceding section since the two end-clocks were 
located on Earth. With our simplifying assumptions, it 
is now a constant differing from unity. Its constancy 
implies that the time delay function may be recovered by 
integrating the Doppler shift with respect to time. 

Explicit expressions for the time delay are obtained in 
appendix ^ For propagation in a static and isotropic 
metric, the time delay only depends on the ratio g r r/goo, 
that is also in the linearized approximation on the grav- 
itation potential <I> . At first order, it is given by IjAlOf) 

rn rf j r 



2 1 *°(r)-^= 
lp - p z 



(59) 



r\2 is the Minkowski spatial distance, i.e. the time delay 
in the absence of gravitational perturbation, 



r\2 



rr 



\Ai 2 + r 2 2 - 2rir 2 cos(( y 9i - tp 2 ) 



(60) 



p is the impact parameter and its time variation is given 
by lHH) 



dp ^Jr x 2 - p 2 ^r 2 2 - p 2 . 

-JT = <Pl 

Otl T\2 



(61) 



(pi is the angular velocity of Earth on its orbit, that is also 
the cause of variation of p in the deflection experiment. 

We now discuss anomalies of the time delay, deflection 
angle or Doppler shift, due to deviations of the potential 
$° from its standard Einstcinian form. To this aim, it 
is convenient to decompose <£>° into a Newtonian part, 
isolating the singularity at r = 0, and a residual part 
expected to become significant at large distances 

r $o (r) ^ _G_M M.2 C o (r) 
c c z 

We consider light rays passing near the solar disk and 
assume a smooth variation of (°(r). Neglecting terms of 
the order or smaller than (p/ri) 2 or (p/^) 2 , we obtain 
in appendix^ the following approximate expression for 
the anomalous time delay 



cST 



G N M 4nr 2 
Sj(p) In — — 



(63) 



2M 



(°(r)rdr+ / (°(r)rdr 



where £7 is the anomalous contribution to an effective 
Eddington parameter 7 



2(G°-G N )-C°(P)P 2 



G 



N 



(64) 



We then deduce the anomalies for the deflection angle 
(5-0 and Doppler signal 5y + 



Sip 
5y+ 



r 2 djcST) 
r\ + r 2 dp 
d(6T) nfa 



dti 



Si/) 



(65) 



Expressions l|63|) and (|65|l contain the same information, 
as 5ip corresponds to the derivative of ST with respect 
to p and 5y+ to its derivative with respect to t\. Note 
that the second term in expression l|63|l is constant and 
does not contribute to 1651) . Both expressions remain 
valid when the emitter is far outside the solar system, 
provided its distance r 2 be replaced by the cutoff distance 
r\ (see the discussions at the end of section Inland in the 
appendix) . 
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The usual PPN expressions are recovered when (°(p) 
vanishes, so that 5j(p) reduces to a constant. The time 
delay and deflection angle thus take simpler forms 



cWppn 

SljjpPN 



G N M 4nr 2 

<i>7o — r~ ln — r~ 

cr p z 

G N M 



>2 



c 2 p 



(66) 



with (5-7o identified as the anomalous part (7 — 1) of the 
Eddington parameter 7 



2(G° - G N ) 



G 



^ (67) 



G 



N 



The main difference with the more general expressions 
(|63|l and (|65|l is that the dependence of £7 on p poten- 
tially present in these expressions has disappeared from 
the PPN expressions l|66|) . This points at a remarkable 
consequence of the extended framework, namely the fact 
that a deflection test may contain more information than 
a constant anomaly S^q. Indeed, a continuous measure- 
ment of the deflection angle during the occultation by the 
Sun may in principle provide one with the dependence 
(|65|l of the deflection angle versus the impact parameter. 

Using again the fact that the anomalous part 5&n of 
the first potential <&jv has to remain very small, this en- 
tails that precise information on the second potential $ p 
could become accessible from a confrontation to obser- 
vations of the formula l|65(l . An observation of a varia- 
tion of £7 with p would constitute a clear confirmation 
of the extended framework proposed in this paper. In 
any case, deflection data have to be confronted against 
Pioneer data since they are determined by the same func- 
tion <f>p. This point is discussed in a more quantitative 
manner in the next section. 



VII. DISCUSSION 

We now collect the results obtained in the previous sec- 
tions to confront available observational or experimental 
data to the extended framework studied in this paper. 
We first write down quantitatively the stringent bounds 
on the anomalous part S^n of $jv which are drawn from 
planetary data. We then focus the attention on the con- 
tribution of the second potential <I>p, considering suc- 
cessively the problems of Pioneer-like anomalous accel- 
erations and Eddington deflection experiments. For the 
sake of making the discussion more precise, we occasion- 
ally use the simple model (|19fl as a benchmark, keeping 
however in mind that the extended framework naturally 
accomodates a more general spatial dependence of the 
gravitation potentials <&at and <&p. 

Constraints on 6<&n may be drawn from a compari- 
son of the orbital frequencies of different planets, pro- 
vided the radius of the orbit is measured through opti- 
cal ranging, that is independently of the measurement of 
the orbital period |57|. This has been done for example 



for Mars, using telemetry data from Viking probes |43| . 
Constraints on the anomaly of the Newton force 9 r $jv 
are then deduced by using the third Kepler law f4*I)l. This 
argument was alread y p resented in section XI-B of [Icj . 
It was developed in |4q . using the most recent bounds 
for violations of Newton law at distances of the order of 
planetary orbital radii (5^ ■ 

We write here the resulting bound as a difference of the 
ratios d r 5$>N /d r (j> at the radius of Mars (tm ~ 1.5AU) 
and Earth (r E = 1AU) 

d r 6<f> N (r M ) d r 6$ N (r E ) < 10 _ 10 ^ 
d r (/)(rM) d r (f)(r E ) ~ 

In the simple model l|19fl . this is read as a bound on the 
anomalous parameter £jv 



\CnM\ < 10- 10 x ~ 5 x 10- 13 ms- 2 



(69) 



This value is small enough to discard any possibility that 
the Pioneer anomaly could be explained from a long- 
range modification of the Newton law 10]. In other 
words, an adhoc modification of first potential $jy de- 
signed to fit the Pioneer anomaly would produce an ef- 
fect on outer planets more than 1000 times too large 
to remain unnoticed |48|. Now, we have shown in the 
present paper that the second potential $p leaves the 
third Kepler law unaffected whereas it has an influence 
on Doppler tracking of probes with eccentric motions. 
This opens free space for a gravitational interpretation 
of the Pioneer anomaly without conflict with planetary 
data |3^ |. 

From now on, we disregard the anomalous part S&n 
of the potential $ n which has been shown to have a neg- 
ligible effect. We focus the attention on the contribution 
of the second potential <E>p to the Pioneer anomalous ac- 
celeration (|56|) which is now read 



So 



drn 



dr 2 
ds 2 



cfr 2 
' ds 2 2 



(70) 



The term proportional to the kinetic energy in l|70|l is 
found to dominate the other one which scales as the grav- 
itational energy of the probe. Keeping only the former, 
we note that the anomalous acceleration is proportional 
to the kinetic energy of the probe, which is a remark- 
able prediction of the framework studied in the present 
paper. Should we have at our disposal data points as- 
sociated with probes with very different kinetic energies, 
this prediction could be used to confirm or infirm the 
framework. 

In the simple model <|19fl . the anomalous acceleration 
(j?0"|l is read as 



8a ~ 2 ( P M 



G P M \ v P 



(71) 



In this equation, the acceleration GpM/rp 2 is much 
smaller than the Newton acceleration GnM/tp 2 since 
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the ratio \Gp/Gn\ ~ |<$7o| is much smaller than unity. 
It is also smaller than the acceleration (pM if the latter 
is to be identified with the Pioneer anomaly (see below). 
As the Pioneer 10/11 probes follow an escape trajectory 
towards the outer solar system with their velocities re- 
maining practically constant 



v P ~ 1.2 x 10 4 ms _1 



Vp 



1.6 x 10" 



(72) 



the anomalous acceleration 5a is found to be approxi- 
mately constant when the distance to the Sun varies, in 
conformity with one of the most striking properties of Pi- 
oneer observations [Tol |. Let us however emphasize that 
this conclusion has been obtained after various simplifi- 
cations recalled in the present section whereas the more 
general expression written above may in principle be con- 
fronted against more detailed data. 

If we now make the jump of identifying the anoma- 
lous acceleration 6a recorded on the Pioneer 10/11 probes 
with the effect (|71|) induced by the second potential <I>p, 
we get an evaluation of the constant C,pM 



8a 



2( P M 



8 x 10" 10 



ms 



( P M ~ 0.25ms" 



(73) 



Note that the positive sign of 5a is associated with a blue 
shift of Doppler frequencies, with the same convention as 
in ^(j- Hence, the identification l|73[l entails that the con- 
stant £p is positive. Comparing with (|69fl . one also notes 
the important difference of magnitude between the linear 
coefficients and C,p in the two gravitation potentials. 

At this point, it is worth coming back to the condi- 
tions of validity of the linearized approach which were 
discussed at the end of section ITTT1 Putting the numerical 
values obtained in l|69[l and l|73[l in these discussions leads 
to an a posteriori justification of the perturbative treat- 
ment performed at first order around Minkowski metric. 
For escape trajectories towards the outer part of the solar 
system, second order terms in the gravitation potentials 
are negligeable when compared with terms proportional 
to linearized potentials and kinetic energies. This stands 
in contrast with the assumption of a comparable magni- 
tude for kinetic and gravitational energies on which the 
PPN framework is built up. While such an assumption 
is well suited for describing the circular motion of plan- 
ets, which justifies its key role in the PPN framework, it 
is also clear that it does not hold for the discussion of 
Pioneer-like eccentric motions. 

We now come to the discussion of deflection experi- 
ments which are determined by the same potential <f>p(r). 
To this aim, we first rewrite the definition of the anoma- 
lous part 5j of the Eddington parameter, using the fact 
that <5$ n may be disregarded. We decompose $p as was 
done for $° in (|62|l . thus separating a part singular at 
r = and a regular part describing the long range effect 



r$p(r) 



G P M M 



r 2 ( P (r) 



(74) 



As (n{t) is negligible, (p{r) is identical to —2(°(r) with 
£°(r) the quantity used in the preceding section. For rays 
passing near the surface of the Sun, we then obtain the 
anomalous time delay l|63l) 



c5T 



. . ,G N M 4nr 2 



M 



C,p{r)rdr + / C,p{r)rdr 



5j(p) = <57o + 



(p(p)p 2 
2G 



N 



(75) 



We have assumed a smooth variation of (p( r ) around 
r = and neglected terms which are not amplified near 
the occultation. We then rewrite the anomaly l|65|l of the 
deflection angle as 



5ip 



r 2 G N M f25j{p) d6-y{p). 4rir 2 

-in — - — 



n + r 2 



V P 



dp 



(76) 



The PPN expressions are recovered when Cp(p) vanishes, 
so that 5j(p) reduces to the constant 570. 

As already discussed, the specific feature of the ex- 
tended framework is the dependence of £7 on p which re- 
veals the long range variation of the potential $p. This 
specific feature can be looked for in deflection experi- 
ments by keeping trace of the whole functional depen- 
dence of the anomalous deflection angle dip during an 
occultation by the Sun 



pdip 5~/{p) 



psi>s 



pd5-/(p) 4nr 2 . . 

-Idp- ln — (7?) 



In order to define manifestly dimensionless quantities, we 
have introduced the Sun radius ps and the standard de- 
flection angle -0s = ^GnM/c 2 ps for a ray grazing the 
Sun. For most experiments, r\ has to be replaced by 
the Sun-Earth distance rp and r 2 by the Sun-probe dis- 
tance. When the emitter is a source far outside the solar 
system, r 2 has to be replaced by the cutoff distance r\. 
Should the precision of the measurement be sufficient, we 
could reconstruct the p-dependence of #7 and, then, the 
r-dependence of Cp and <I>p. 

We want to stress that the purpose of extracting the 
p-dependence of the effective parameter ^7 is different 
from that of existing analysis which aim at obtaining a 
global value of £7 on the whole range of available values 
of p. Formula lj77|l suggests that deviations should tend 
to vary with the impact parameter, being larger when 
the impact parameter increases, that is also when the 
deflection angle decreases, so that the detection strat- 
egy cannot be the same as that used in existing analysis. 
Consequently, the constraints drawn for example from 
Cassini j3^| or VLBI measurements [3^ cannot easily be 
translated into a constraint on the p-dependent contri- 
bution to 5^. We can only make a guess whether or not 
these measurements can produce information of interest 
for our purpose by proceeding as follows. 
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If we consider the simplest model (|19ll with (p constant 
over the whole range of distances from the radius of the 
Sun to the size of the solar system, we may use the value 
(|73|l of (p deduced from Pioneer data as a first guess for 
the Cassini experiment. We thus obtain from (|77ll 

pStp 5^o (pMp 2 4rir 2 

— r - 2 — r ln — 2~ ( 78 ) 

p s ips 2 c z psws e.p 

This quantity varies by a few 10 -3 when p goes from 5 to 
10 p $, which should probably be visible in Cassini data 
[38j . if looked for. Emphasizing once again that the sig- 
nature we are interested in is a potential variation of 5"f 
with p rather than a mean value of £7, this result cer- 
tainly pleads for a reanalysis of Cassini data |38j aimed at 
the observation of such a variation. Such an observation 
would constitute a clear indication for a non null value 
of Cp, to be then confronted to the value deduced from 
Pioneer data. Different values could also be obtained in 
the two experiments since (p is in fact a function of r 
which could vary when r goes from the radius of the Sun 
ps ~ 0.7 x 10 9 m to the distance explored by the Pioneer 
probes rp ~ 1.2 x 10 13 m. 

VIII. CONCLUSION 

We have shown that Einstein gravitation theory pos- 
sesses natural extensions characterized by two running 
gravitation couplings replacing the single Newton con- 
stant or, equivalently, two gravitational potentials $at 
and $p replacing the usual Newton potential </>. The 
running couplings may depend on the length scale and 
differ in the two sectors of traceless and traced tensors, 
which have different conformal weights. 

The good agreement between planetary data and Gen- 
eral Relativity essentially means that the long range mod- 
ification of the first potential $at has a negligible ef- 
fect. This conclusion has been drawn from the discus- 
sion of third Kepler law for circular orbits. It would 
also be worth discussing the perihelion precession but 
this study, which involves nonlinearities of the gravita- 
tion theory, is postponed to a forthcoming analysis. The 
linearized treatment performed in the present paper en- 
compasses the phenomena for which the kinetic energy 
is much larger than the gravitational one, that is also 
those usually ascribed to the 7 parameter in the PPN 
formalism. This includes the case of spacecrafts follow- 
ing escape trajectories towards the outer solar system as 
well as deflection experiments on light rays. 

The second potential $p has been shown to have a 
crucial effect on Doppler tracking of probes with highly 
eccentric motions, allowing one to propose a gravitational 
interpretation of the Pioneer anomaly without raising a 
conflict with planetary tests. This interpretation comes 
with the crucial prediction of a magnitude of the anoma- 
lous acceleration proportional to the kinetic energy of 
the probe. This prediction cannot be confronted against 
the data available today for the Pioneer 10/11 probes 



which were endowed with very close kinetic energies, but 
it could be the subject of dedicated analysis of Doppler 
tracking data stored during their flight to the outer solar 
system. In particular, a study of the fly-by sequences of 
the probes might provide interesting clues. 

While the presence of <&p at scales of the order of the 
solar system size is tested by Pioneer-like Doppler track- 
ing, its variation in the vicinity of the solar radius can be 
studied through deflection experiments. We have shown 
that $p may thus be seen as promoting the anomalous 
Eddington parameter (7 — 1) to the status of a function 
of the impact parameter p. Evidence of a variation of 7 
would provide a direct validation of the extended frame- 
work. This aim might be attainable through a reanalysis 
of exis ting data, in particular those of the Cassini exper- 
iment 38]. On a longer term, it can also be reached by 
high accuracy Eddington tests (see for example the LA- 
TOR project [5!|) or global mapping of deflection over 
the sky (see for example the GAIA project 60] ). The 
former would lead to an accuracy on the measurement 
of 7 certainly sufficient to see the effect if its order of 
magnitude is given by the numbers deduced from the Pi- 
oneer data in the preceding section. Meanwhile the latter 
would allow one to reconstruct the dependence of 7 ver- 
sus the impact parameter p and then, the dependence of 
<&p versus r. Should these future projects conclude to 
the absence of such a dependence, this would anyway re- 
sult in improved constraints on a family of natural metric 
extensions of General Relativity, confirming its validity 
for describing gravity in the solar system. 



APPENDIX A: TIME DELAY FUNCTION 

In this appendix, we determine explicit expressions 
for the time delay function T which characterizes light 
propagation. T is a two-point function derived from the 
phase, that is also, in the semi-classical limit, the action 
S which satisfies Hamilton- Jacobi equation. 

Stationarity and isotropy of the metric entail the con- 
servation laws H32|) . It follows that the action S can be 
written in terms of a phase shift S r which is a function 
of r only (see fti'6\). The phase shift obeys an Hamilton- 
Jacobi equation i|30fl rewritten here for light propagation 
(massless field) 

{drSrf = ~ ^ (Al) 

3oo c z r 1 

It can be given the following integral form determined 
by the conformally invariant ratio g rr /9oo °f t ne metric 
components 

Sr = ±J f drJ~^- 2 9 —-- 2 (A2) 

J v J c goa r 

The same solution is equivalently written as a character- 
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istic equation for the associated null geodesies 
±l/r 2 



dip 
dr 



E 2 g rr _ J_ 
J 2 c 2 g 00 r 2 



djd r S r (A3) 



A null geodesic is essentially characterized by the radial 
coordinate p at closest approach of the Sun. Equation 
(|A2(I is then solved in terms of the reduced function 

w{x) = — {p) — (r) , x = - (A4) 
grr goo r 

Fixing by convention the value of S r to be zero at closest 
approach, its explicit expression is obtained as 

S r {r) = -Js w (-J , ip < ip 
S r (r) = +Js w , ip > (p 



,(x) = [ ^y / w(x) - X 2 



(A5) 



The geometry of the null geodesic is then determined by 
writing boundary conditions on its enpoints 



<f{r) 



V r 



Vo, 



P < Po 



p(r) - Cw y-J = (fo, 
c w (x) ee -dj (Js w (x)) 



P>P>0 



dx 



x y / w(x) - X 2 



P2 + c lu \ ~j =(p = (px-c w (~^J (A6) 

Note that equations (|A6|I simultaneously determine the 
distance of closest approach and the angles ipo as func- 
tions of the endpoints. 

The time delay, that is the elapsed time during propa- 
pagation of a constant phase from (ri, ipi) to (r2, ^2), is 
then deduced from the phase shift (|A5|I 



T12 = T(rx,(pi,r2,(p7;) 

= TT {J ('Pi - V2) + S r (n) - SV(72)) 

hi 



(A7) 



grr , sP IP 

LP)- ['Pi - + Sw — I + s. 

goo c \ \nj \r 2 



The wavevector, tangent to the null geodesic, can be ob- 
tained from the endpoint variations of the time delay 



di(cT 12 ) = - 



goo (n) 



(A8) 



9rr(n) 

x (cos(-0i — tpi)dr\ + sin(^i — ipi)r\dipi) 

Here d\ represents a variation of the endpoint (r\,ipi) 
with the other one kept fixed and ip\ denotes the line 
of sight angle, that is the angle associated with the light 
propagation direction. Variations with respect to the sec- 
ond endpoint may be recovered by symmetry. Equations 



written up to now hold at any order in the gravitational 
perturbation. 

In the absence of gravitational fields, cl\i is given by 
the spatial distance 7*12 evaluated in Minkowski metric 
(see eci. (|60[l '). Meanwhile null geodesies are straight lines 
and the line of sight angle is preserved under propagation 



1p2 =1pl= 1p(Q) 



(A9) 



We now simplify the expression l|A7(l of the time delay 
function up to first order in the potential $° 



\J 'r 2 — p 2 



- 2 



(A10) 
rdr 



y/ 'r 2 — p 2 



Expanding relation l|A8|l up to first order in $°, one ob- 
tains the deflection angle Ai/j at the same order 

= Vi-V'(o) (All) 
rfi(cT 12 ) = (l-2$ (r 1 ))dr 12 +r 12 A^ (0) 

Note that the derivative of the time delay along the line 
of sight is simply related to the gravitation potential <I> . 
Equation (|A11|) involves the unperturbed angle "0(0)? or 
equivalently the impact parameter p, at lowest order 



dip 



-\Jr 2 2 - P 2 dip {0) 

vVl 2 - p 2 {dxp( ) - dipt) + ydr 



(A12) 



This entails that the deflection angle is simply related 
to the derivative of the time delay with respect to the 
impact parameter, at constant radial coordinate r\ 



Atp = - 



y/r 2 2 - P 2 Q(cT) 
r\i dp 



-(i-2$°(n)) 



yVi 2 - p 2 



(A13) 



Expression 1A10I) shows amplification of the time de- 
lay near occultation. In order to isolate the singularity 
at vanishing impact parameter we decompose the poten- 
tial < I >0 into its Newtonian singular part and a residual 
representing the long distance modification, according to 
(|62|l . We then rewrite it as 



+ r 2 (c°M- C °(p)) 



M 



(A14) 



In order to allow for a simple comparison with the PPN 
framework, we now introduce a function "/(p) which gen- 
eralizes the Eddington parameter (see ea. l|t)4|0 



l(p) 



2G°-CW 
Gn 



(A15) 
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The time delay IjAlOf) then takes the explicit form 

cT 12 = ru (A16) 

- (i + t(p)) — — ln 



cr ri + r i2 + r 2 



M 
72" 



((2C°(r 1 )-C°(p))r 1V / ^7 
+ (2C°(r 2 )-C°(p)) r 2 V^~?) 
+ ^ (f V^?d (r(C°(r) - C°(p))) 

+ y/*=?d(r(<?(r) -C<V)))) 

The logarithmic part has the same form as in the PPN 
formalism, but with 7 now a function of the impact pa- 
rameter p. This dependence is determined by the long 
range behavior of the potential <£>°, that is also by the 
function £°. Assuming that £° is a regular function, the 
last terms of equation (|A16(1 may be expanded in pjr\ 
and p/r2, so that the time delay (jAlOfl may be rewritten, 
up to terms equal to or smaller than {p/ri) 2 or (/o/r 2 ) 2 



c7i2 = r%2 



(A17) 



(1 + 7(P)) 



-In 



r 1 + ria + r 2 



In the same context, the deflection angle (|A13|) may be 



rewritten up to terms of order p/r\ or p/r 2 
r 2 d(cT) 



ri2 dp 

r 2 G N M d 



(A18) 



n + r 2 



dp 



(l+ 7 (p))ln^p 



The logarithmic divergence with respect to r 2 in equa- 
tion l|A18(l signals a failure of the perturbation expansion 
leading to equations l|A10|l and IjAlip . due to the fact 
that the gravitational potential cannot be considered as 
weak over large distances. To fix this problem, one must 
come back to exact expressions given previously for the 
deflection angle. As discussed in section (HIJ, the linear 
behaviour of the long range contribution to <J>° can only 
be an approximation within the solar system of the true 
potential $° which should remain weak over much larger 
distances. This behavior is thus limited to a maximal dis- 
tance ta, beyond which the potential recovers a purely 
Newtonian form. For instance, the cut-off may take the 
form of a mass parameter in a Yukawa potential with a 
range ta larger than the size of the solar system . The 
effect of the anomalous part of <J>° from this range to 
very remote light sources may then be neglected, and the 
deflection angle obtained by replacing r 2 by r\ in i|A18(l 



G N M d 
c 2 dp 



(l + 7 (p))ln 



4nr A 



(A19) 
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